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Abstract 


The usage of Computational Fluid Dynamics (CFD) in noise prediction typically has been a two part 
process: accurately predicting the flow conditions in the near-field and then propagating the noise from the 
near-field to the observer. Due to the increase in computing power and the cost benefit when weighed against 
wind tunnel testing, the usage of CFD to estimate the local flow field of complex geometrical structures 
has become more routine. Recently, the Benchmark problems in Airframe Noise Computation (BANC) 
workshops have provided a community focus on accurately simulating the local flow field near the body 
with various CFD approaches. However, to date, little effort has been given into assessing the impact of 
the propagation phase of noise prediction. This paper includes results from the BANC-III workshop which 
explores variability in the propagation phase of CFD-based noise prediction. This includes two test cases: 
an analytical solution of a quadrupole source near a sphere and a computational solution around a nose 
landing gear. Agreement between three codes was very good for the analytic test case, but CFD-based noise 
predictions indicate that the propagation phase can introduce 3dB or more of variability in noise predictions. 

Nomenclature 


English: 


Greek: 

c 

speed of sound 

S(x) 

Dirac delta function 

f 

surface of a body (/ = 0) 

hj 

Kronecker delta 

Ft 

compact line dipole source term 

P 

fluid density 

G 

Green’s function 

a 

standard deviation 

H(x) 

Heaviside function 

T 

retarded time 

Li 

surface dipole source term 

Subscript: 

M 

Mach number 

/ 

incident field 

7li 

outward directed unit normal vector 

n 

unit normal direction 

p' 

acoustic pressure 

S 

scattered field 

Pij 

compressive stress tensor 

oo 

freestream quantity 

Q 

monopole source term 

Superscript: 

r 

radiation vector 

/ 

perturbation quantity 

S 

integration surface 

- 

generalized derivative 

t 

observer time 

Symbol: 

T- 

J y 

quadrupole source term 

□ 2 

wave operator, 1 /c 2 OG d ‘ 2 / 

Ui,Vi 

flow and surface velocity 
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I. Introduction 


The study of noise from aircraft has historically progressed on three fronts: analytical theory, wind 
tunnel or field measurements, and computer-based applications including Computational Fluid Dynamics 
(CFD). Typically, utilizing CFD to predict noise involves a large calculation that captures the near-field 
fluctuations including acoustics and a post-processing algorithm that propagates the acoustics from the 
near-field to observer positions. Separating the computation into two parts allows the propagation phase of 
noise prediction to use reduced order models in order to save time and computer resources. However, there 
are several different approaches to noise propagation, each with benefits and assumptions that affect the 
overall result. This leads to uncertainty and variability that can be directly associated with the propagation 
phase of noise prediction. 

The propagation of noise from a sound generating region to an observer position involves many problem 
areas such as scattering obstacles and medium heterogeneities such as shear layers. There are several different 
approaches to the propagation phase of noise prediction including Integral Methods (such as Kirchhoff 
extrapolations 1,2 and Acoustic Analogies 3,4 ), Ray Methods (RM), 5 Boundary Element Methods (BEM), 6,7 
Equivalent Source Methods (ESM), 8,9 and Computational Aero Acoustics (CAA). 10 While there has been 
effort in each area to categorize and validate approaches, there has not been a concentrated effort to compare 
the range of applicability and trade-offs across the methods. 

Recently, the Benchmark problems for Airframe Noise Computations (BANC) workshop has focused on 
the validation of CFD-based approaches for noise prediction of airframe sources. 11 The BANC workshop 
includes eight categories: the first seven focus on near-field computations of the flow field around airframe 
structures, such as the landing gear, flap, and trailing edges; Category 8 focuses on assessing the variability in 
the propagation phase of CFD-based noise prediction. This paper presents results from BANC-III Category 8. 
This includes applying different propagation approaches to two subcategories: the scattering of a quadrupole 
source by a sphere and the propagation of noise from a simplified nose landing gear. A summary of the BANC- 
III Category 8 participant’s propagation methods is included, after which the result from each subcategory 
is shown. Finally, a Conclusion and Future Work section is presented that includes proposals for additional 
sub categories for follow-on efforts. 


II. Propagation Methods 

The noise propagation portion of CFD-based noise prediction involves using CFD to predict the time- 
dependent near-field hydrodynamics and acoustics and a propagation algorithm to calculate a pressure time 
history at an observer position. Outside the near-field where the noise is generated, there can exist significant 
influences on the noise propagation, such as shielding/scattering bodies or flow heterogeneities. These reside 
in a mid-field region. Once outside the mid-field region, noise typically propagates through a homogeneous 
medium to an observer. Figure 1 shows a comparison between a simple propagation problem and a more 
realistic one. Figure la shows propagation with no body interactions and an immediate transition from 
near- to far-held. Figure lb shows propagation with potential shielding and scattering interactions, as well 
as heterogeneous propagation effects in the mid-field. There are several different approaches to calculating 
the noise propagation, summarized in Figure 2. Each of the approaches must be provided a solution to 
the near-held region, including hydrodynamics and acoustics, and can provide a prediction at an observer 
position in the near- or far-held. Once a pressure time history at the observer is determined, post processing 
algorithms are applied so that a direct comparison to wind tunnel measurements can be performed. Due to 
participation in Category 8 of the BANC-III workshop, for brevity, only three approaches are included in 
this study: Integral Methods, Boundary Element Methods, and Computation Aero Acoustics. Following is a 
brief explanation of each. 
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Figure 1. Summary of regions of CFD-based noise prediction, a) Noise generating region placed in isolation 
where no bodies influence noise propagation, b) Body exists near noise generating region influencing the noise 
propagation. 


A. Integral Methods (IM) 


Over the past 30 years, Integral Methods (IMs) 
have provided a relatively computationally inex- 
pensive method for propagating noise from a near- 
field region to an observer location. Historically, 
the most popular IMs are based on a reformula- 
tion of the Navier-Stokes equation using Generalized 
Functions, 12 called the Ffowcs Williams and Hawk- 
ings (FW-H) equation. The FW-H equation utilizes 
Light hill’s acoustic analogy to formulate an inho- 
mogenous wave equation where source terms are de- 
fined as discontinuous surfaces and a volume region 
external to the surface. The FW-H equation, Eq. 1, 
includes two surface source terms and one volume 
source term. The surface source terms are identified 
by the delta function, S(x), and the volume source 
term is identified by the Heaviside function, H{x). 

The source terms on the right hand side of the FW- 

H equation often are called the monopole, dipole, and quadrupole terms respectively. The surface /, shown 
in Fig. 3, can be impermeable, such as a rotor blade surface, or permeable, such as a computational surface 
surrounding the entire rotor which provides flowfield information. The general source terms are shown in 
Eq. 2, and the simplified source terms valid for impermeable surfaces are shown in Eq 3. 



Figure 2. Summary of coupled CFD and propagation 
components of CFD-based noise prediction. 
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Integral methods often ignore the quadrupole term with the assumption that all noise sources and non- 
uniform flow propagation effects are included within and thereby captured by the surface terms. IMs also 
often apply a free-space Green’s function by assuming a uniform propagation region; therfore, IMs cannot 
account for the influence of bodies or medium heterogeneities. These assumptions often lead to inaccuracies in 
the prediction. However, usage of integral methods have provided a very reliable prediction of the propagation 
with a relatively small amount of computational effort. 
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Figure 3. Description of a surface in motion by /(x, t) = 0, n = /(x, t) < 0 inside the surface, and /(x, t) > 0 

outside the surface. 


B. Boundary Element Methods (BEM) 

Boundary Element Methods (BEMs) discretize the surface of any shielding/scattering body and combine the 
convective wave equation, shown in Eq. 4, with a Green’s function, shown in Eq. 5. A useful equation for 
the pressure can be obtained by subtracting the pressure times Eq. 5 from the Green’s function multiplied 
by Eq. 4. If the pressure is assumed to be comprised of an incident field and a scattered field, then Eqs. 6, 7, 
and 8 provide integral relations for the two components of the pressure. 7 The incident field is calculated 
similar to an IM, by assuming a free-space Green’s function. The scattered field is calculated by integrating 
surface elements along the scattering body, where each element contributes to the scattered field accounting 
for its own incident field from all other elements. In this way, BEMs can account for a shielding/scattering 
body; however, BEMs have yet to include effects such as flow heterogeneities. Similar to ESMs, BEMs have 
had difficulty at higher frequencies due to a need for reduced element and time step size, but as computer 
power increases BEMs have gained significant traction. 
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C. Computational AeroAcoustics (CAA) 

Computational AeroAcoustic (CAA) methods rely on a discretized form of the linearized Euler equations 
to account for both acoustic reflection/diffraction by impermeable obstacles and refraction by flow hetero- 
geneities. This is in contrast to other techniques that can only model obstacles (such as the BEM or ESM), 
or far-held propagation (such as an IM). A typical situation where an acoustic hybrid approach can be 
advantageously applied is external noise problems where the propagation phase occurs within a complex 
environment, such as the airframe noise emission by aircraft appendages (e.g. landing gear, etc.). Indeed, 
once their generation has been properly simulated, acoustic waves may then be transferred to a CAA solver, 
for the latter to propagate them through the mid-field, while accounting for all the installation effects in- 
duced by either the airframe elements (e.g. rehection/diffraction) or the airflow that surrounds the latter 
(e.g. convection/refraction). 
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CAA methods usually rely on an arbitrary splitting of the usual physical quantities (primitive variables) 
into both a perturbed (subscript p) and a mean (or background, subscript o) component. 

q = q c + % (9) 

where each q quantity is composed of its associated density, p, velocity, v, and pressure, p. Invoking 
the classical small perturbation hypothesis, the background field, q 0 = (p G , v G , p 0 ), can be considered as 
insensitive to the evolution of its perturbed counterpart, q p = ( p p , v p , p ), as long as the latter remains 
sufficiently small in terms of amplitude. This allows the derivation of an equation governing solely the 
perturbed field, q , whilst accounting for its background counterpart, q G : 

<9tUj,(q 0 , q p ) + VF p (q 0 , q p ) = 0 (10) 

In the above equations, u p stands for the perturbed flow-field vector, and F p stands for its associated 
convective flux matrices. For conciseness, these terms are not detailed here but can be found in Ref. 10. 
Equation 10 is then discretized in space and time and solved numerically; the u p term being advanced in time 
while accounting for the mean flow and other effects (boundary conditions, etc). Note that the space/time 
discretization is usually handled with the help of high-order schemes, optimized for minimal dissipation and 
dispersion of the numerically propagated acoustic waves. 

III. BANC-III Category 8 

The goal of Category 8 of the BANC workshop is to quantify variability in the propagation phase of 
CFD-based noise prediction; this includes resources required, computation time, and accuracy of the result. 
Toward that end, two subcategories were created to quantify the variability. Subcategory 1, Scattering of 
Noise by a Sphere, was created as a simple test case to quantify the execution time, resources, and accuracy 
of each approach compared to an analytical solution. Subcategory 2, Nose Landing Gear, was created to 
assess the variability in the results when provided a common CFD solution of a complex noise generating 
region. 

Participants 

There were four participants in the BANC-III Category 8 effort: JAXA’s upacsAcoustics IM solver, 13 NASA’s 
ANOPP2-F1A IM solver, 14 ODU’s TD-FAST BEM solver, 15 and ONERA’s sAbrinA CAA solver. 7 Following 
is a brief description of each code. 

1. upacsAcoustics 

The JAXA’s upacsAcoustics code is an Integral Method solver that utilizes a convective form of Farassat’s 
Formulation 1A to solve for the acoustic pressure in double precision at an observer location provided the 
local flow conditions around a body in uniform, rectilinear flight. Only the monopole and dipole terms are 
included. The quadrupole term is ignored with the assumption that the freestream Mach number is low. 
The code uses a second-order spatial and temporal central differencing scheme with a first-order retarded 
time calculation. 

2. ANOPP2-F1A 

ANOPP2-F1A is the Formulation 1A solver incorporated in the Aircraft NOise Prediction Program 2 
(ANOPP2). 14 ANOPP2 is a mixed- fidelity acoustic framework, capable of combining high-resolution meth- 
ods, such as those based on CFD computations, and low-resolution empirical methods, such as those based 
on experimental data and historical methods, for the application of acoustic research and Multi-disciplinary 
Design Analysis and Optimization (MDAO). 16 ANOPP2-F1A utilizes a second-order central differencing 
scheme for spatial and temporal derivatives and a first-order retarded time calculation, both in single preci- 
sion. ANOPP2-F1A is capable of predicting the near- and far- field acoustic pressure at an observer location 
provided the appropriate flow data on an input acoustic data surface which may be undergoing complex 
motions. The acoustic data surface may be node- or cell-centered, rigid or deforming, 17 structured or un- 
structured and may also be undergoing complex motions such as those found in rotorcraft applications. In 
addition to Formulation 1A, ANOPP2 also includes Formulations l, 18 2B, 19 Gl, G1A, 20 and VIA. 21 
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3. TD-FAST 


TD-FAST (Time-Domain Fast Acoustic Scattering Toolkit) is a time domain boundary integral equation 
solver. 7, 15 It has the distinct advantage that acoustic scattering and shielding solutions at all frequencies 
within the numerical resolution can be captured by one single time domain computation in double precision. 
In addition, broadband sources and time dependent transient signals can be simulated and studied directly 
in TD-FAST. It also can be used to couple with nonlinear CFD simulations where many frequencies are 
generated. Solution of the boundary integral equation in TD-FAST requires only unstructured surface 
triangular and quadrilateral elements. Computational complexity in TD-FAST is reduced by the Time 
Domain Propagation and Distribution (TDPD) algorithm that has been developed based on the delay- and 
amplitude-compensated acoustic field. Computation in TD-FAST is further accelerated by a use of general 
purpose GPU computing. TD-FAST has been implemented in a multi-CPU and multi-GPU platform. 

4- sAbrinA 

sAbrinA 10,22 26 is a structured, time-accurate CAA code that solves either the full or the linearized Euler 
equations in a conservative and perturbed form (with a splitting of the complete variables into a ‘frozen’ 
mean flow and a ‘fluctuating’ perturbation). The solver employs high-order, finite- difference (FD) operators, 
as well as a third-order, multi-stage, Runge-Kutta time-marching scheme. The code deals with multi- 
block structured grids with one-to-one interfaces. Finally, the solver includes the usual boundary conditions 
(reflection by solid walls, non-reflecting and free-held radiation, 10,22 etc.), as well some unique to specific 
applications (such as the Non Reflective Interface technique 27 that was extensively used in the present effort, 
for properly forcing the CAA stage with the near-held datasets) . The present effort benehts from a new class 
of hnite difference (FD) schemes that were recently developed. 28 The higher accuracy of these Intrinsically 
Optimized FD (IOFD) schemes produce accurate solutions on courser CAA grids than can be used with 
standard, order-property operators of the same stencil width. 

Subcategory 1: Scattering of Noise by a Sphere 

Subcategory 1: Scattering of Noise by a Sphere consists of a point quadrupole located near a sphere in a 
quiescent, homogeneous medium. The purpose of Subcategory 1 is primarily for use as a code verihcation 
exercise demonstrating the properties and order of the codes. For this subcategory, participants were provided 
a Fortran code that could calculate the acoustic solution on a box around the point quadrupole and sphere, 
including any scattering/rehection caused by the sphere. The solution on this box was then provided to 
the propagation algorithms which were required to calculate the noise on an arc outside the box. The noise 
prediction was then compared to the analytical solution. Participants were asked to provide information 
about the computational resources including architecture, details on processing unit(s), wall-clock time, and 
memory required. Figure 4 shows a schematic of Subcategory 1 including the point quadrupole, the sphere 
which acts as a scattering body, and the six-sided computational surface. 

Subcategory 1 involved three participants: upacs Acoustics, TD-FAST, and sAbrinA. Figure 5 shows a 
snapshot of acoustic pressure predicted by upacsAcoustics on a sphere surrounding the point quadrupole. 
TD-FAST and sAbrinA provided similar data. Table 1 tabulates the different approaches to subcategory 1, 
including the number of observers where a prediction was made, the architecture of the computers utilized, 
and the total amount of memory required to perform the prediction. upacsAcoustics and TD-FAST predicted 
the noise on a similar number of observers, but since sAbrinA is a CAA method, the solution at every cell used 
in the computation is available, thus the number of observers for sAbrinA is significantly higher. The TD- 
FAST and sAbrinA codes required significantly more memory than upacsAcoustics and more computational 
effort. Table 2 summarizes the wall clock time and error of the computations. All three methods produced 
similar error; however, the wall clock time between the methods is very different. Since sAbrinA is a CAA 
method that must calculate the solution on a mesh encompassing the entire propagation region, there are 
significantly more observer positions. The wall clock time per observer per time step is the smallest. TD- 
FAST requires the most computational resources per observer, but this metric is biased because of the 
number of observers associated with the CAA method. In most instances, the solution is not needed at all 
the points in the CAA mesh. 
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Analytical solution at t = 0 


Figure 4. Schematic of Subcategory 1 which includes a point quadrupole near a sphere of radius R and the 
computational surfaces surrounding both. Also shown is the acoustic solution on the computation surfaces at 


t = 0. 



Approach 

^ Observers 

Architecture 

7 ^ Cores 

Memory (Mbytes) 

upacsAcoustics 

IM 

342 

Not Available 

1 

35 

TD-FAST 

BEM 

342 

Nvidia M2090 GPU 

24 

800 

sAbrinA 

CAA 

520000 

Intel 17-3610 QM CPU 

1 

9,182 


Table 1. Summary of participants’ setups in subcategory 1. 


Comparisons of efficiency and speed between im- 
plementations of methods that are significantly dif- 
ferent in approach is difficult. Due to time limi- 
tations, variations in computational setup, such as 
time step size and mesh size were not included in 
this work. In the future, more exploration into the 
potential variation due to these parameters will be 
explored. Even so, the above represents a first step 
in this comparison. In future workshops, additional 
participation will result in a better understanding of 
the variation in computational requirements. Also, 
improvements will be made to Subcategory 1, in- 
cluding a more complex noise calculation (such as a 
scattering body outside the source region) and some 
metric for benchmarking the computational require- 
ments of the predictions. These efforts would result 
in a more consistent and representative metric for 
all methods. 

Subcategory 2: Nose Landing Gear 


z 



Figure 5. Snapshot of instantaneous pressure field pre- 
dicted by upacsAcoustics on sphere of 5R surrounding 
point quadrupole. 


Subcategory 2 was created to quantify the variabil- 
ity in predicted spectra from propagation algorithms 
given a more realistic input. Participants were pro- 
vided a CFD solution of the flow around a nose landing gear geometry (the Airbus-ONERA LAGOON 
configuration used in BANC-II Category 5). 29 31 The solution, provided as an internet download of di- 
mensional flow quantities, included on-body (impermeable) and off-body (permeable) surfaces, the latter of 
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Wall Clock / Observer / Time Step (/i sec) 

Mean-RMS Error 

upacsAcoustics 

20.3 

4.3369e -6 

TD-FAST 

1950 

2.7160e -6 

sAbrinA 

8.6 

1.782e -7 


Table 2. Summary of wall clock time and mean-rms error for participants of subcategory 1. 


which is designed as a nearly cubic interface encompassing the wheels and is composed of 6 panels shown in 
Fig 6. A Fortran 95 code that reads the datasets was also provided to the participants for ease in handling 
the data. The impermeable and permeable datasets contained 20100 and 14390 time steps, respectively, 
both with a common time step spacing of 5e-6 seconds. Participants were permitted to downsample the 
CFD solution in time and/or space as needed to help characterize the dependence of their predictions on 
spatial and/or temporal resolution. Participants were asked to provide the acoustic pressure at 5 far-held 
observer locations (designated FI through F5), 4 near-held observer locations (designed N1 through N4), 
and at an observer at the midpoint (origin) between the wheels inside the surfaces (designated as the origin). 
Supplied results were to include pressure time signals and power spectral densities. The observer locations, 
with the exception of the origin, coincided with measurement locations used in ONERA’s open-jet anechoic 
windtunnel CEPRA19 in Saclay, France. 30 

Subcategory 2 included three participants: two 
IM solvers (ANOPP2-F1A and upacs Acoustics) and 
one CAA solver (sAbrina). ANOPP2-F1A and up- 
acsAcoustics utilized the entire CFD solution (not 
downsampling in space and time); however, sAbrina 
utilized a solution that was downsampled in space 
(1:2 ratio) and time (sufficient to resolve a maxi- 
mum frequency of 3kHz). Figure 7 shows a snapshot 
of predicted acoustic pressure time histories com- 
puted using the permeable surface data. Although 
there are signihcant variations in the three predic- 
tions, they show similar trends. Some of the varia- 
tion can be attributed to the downsampled solution 
provided to the sAbrina code. It is also interest- 
ing to note that even though upacsAcoustics and 
ANOPP2-F1A utilize the exact same formulation 
for noise propagation, there are differences in the 
predictions that are non-trivial. All three methods 
included different observer time ranges and averag- 
ing techniques that complicated comparisons of the 
supplied spectra. Therefore, the spectra were com- 
puted in a consistent fashion based on the supplied 
time history data. A common temporal sample of 
0.069 seconds for predictions based on the imperme- 
able surfaces and 0.042 seconds for predictions based on the permeable surfaces was chosen. This resulted 
in 1666 samples from upacsAcoustics and sAbrinA and 3400 from ANOPP2-F1A. A Hanning window was 
applied to all time histories before calculating a power spectral density. Figures 8 and 9 show the power 
spectral densities from the permeable surface calculations using the full data record. The results are pre- 
sented in a narrow 23.8 Hz bin width. It is difficult to assess the three results using the narrow bin width 
due to the high variability. A longer time history providing sufficient averaging, would result in a smoother 
power spectral density curve; however, due to the limited duration of source time history in the CFD input, 
a longer predicted time range is not possible. Therefore, an FFT averaging scheme was employed where the 
time histories were separated into 5 FFT averaging blocks for predictions based on the permeable dataset 
and 8 FFT averaging blocks for predictions based on the impermeable dataset. Both averaging techniques 
employed 0% overlap resulting in a 120 Hz bin width. 

Figure 10 shows predictions at FI through F5 from the three participants using a permeable surface solu- 



Figure 6. Permeable interfaces extracted from CFD 
solution of flow around LAGooN geometry. 
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Figure 7. Snapshot of acoustic pressure time history predictions at FI from ANOPP2-F1A, upacs Acoustics, 
and sAbrinA codes using permeable surfaces in shape of a box around landing gear. 






Figure 8. Power spectral densities from ANOPP2-F1A, upacs Acoustics, and sAbrinA codes using permeable 
surfaces at microphone locations FI through F5. No FFT averaging was applied. 
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Figure 9. Power spectral densities from ANOPP2-F1A, upacs Acoustics, and sAbrinA codes using permeable 
surfaces at microphone locations N1 through N5 and the origin. No FFT averaging was applied. 


tion as input. For all observer positions, all predictions agree very well with the experiment and each other. 
It is interesting to note that all predictions tend to overpredict around and below 1kHz and underpredict 
in the high frequencies. Since all methods are behaving similarly, this trend indicates a dependence on the 
simulation data used as input rather than some characteristic of the codes. Undoubtedly, the diffusive nature 
of the CFD simulation results in an artificial decay of high-frequency waves (above 2kHz). The sAbrinA 
results decay even faster at high frequencies, most likely because of subsampling the input data, but possibly 
because of diffusion in the calculation itself. Figure 11 shows the predictions compared to measurements 
at the four near-field microphone locations, N1 through N4. The agreement with experiment for all three 
solvers is very good, particularly at the upstream position Nl. For the near-field microphones Nl, N2, and 
N3, all predictions agree very well with the experiment in the low- and mid-frequencies. At the downstream 
location, N4, all codes significantly overpredict the lower frequencies. Figure 11 also shows the prediction 
at the origin from the two IM solvers, ANOPP2-F1A and upacsAcoustics. Predictions using the FW-H 
equation theoretically should predict exactly zero inside the surface and both codes do not, with ANOPP2- 
F1A predicting significantly higher levels than upacsAcoustics. Possible causes for the discrepancy include: 
neglecting the viscous terms in the Lighthill stress tensor, neglecting the quadrupole term, and numerical 
error. The exact cause is unknown, but will be a focus of future efforts. 

Figures 12 and 13 show the far-held and near-held predictions from ANOPP2-F1A and upacsAcoustics 
using the impermeable dataset. For all microphones the high frequencies are underpredicted by all the 
codes, especially at the near-held microphones. For all microphones, far-held and near-held, the low- and 
mid- frequency levels are predicted extremely well by both codes (the predictions by ANOPP2-F1A are 
hidden by the predictions by upacsAcoustics), including the N4 position where the predictions using the 
permeable dataset overpredict the levels compared to the experiment. This suggests that the overprediction 
in the downstream location and in the far-held microphones at the low frequencies using the permeable 
datasets is due to the conditioning of the rear permeable surface and that using permeable datasets in 
the future requires a correction in order to accurately predict in the low frequency. Figure 13 also shows 
prediction from both codes at the origin, which is placed at the center of the hub. Both ANOPP2-F1A and 
upacsAcoustics predict signihcantly high values due to the closeness of the surface when theoretically they 
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Figure 10. Power spectral densities at microphone locations FI through F5 from ANOPP2-F1A, upacs Acous- 
tics, and sAbrinA codes using permeable surfaces. Five FFT averages were applied resulting in a 120 Hz bin 
width. 
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Figure 11. Power spectral densities at microphone locations N1 through N5 and the origin from ANOPP2- 
F1A, upacsAcoustics, and sAbrinA codes using permeable surfaces. Five FFT averages were applied resulting 
in a 120 Hz bin width. 
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Figure 12. Power spectral densities at microphone locations FI through F5 from ANOPP2-F1A and upac- 
sAcoustics codes using impermeable surfaces. Nine FFT averages were applied resulting in a 120 Hz bin 
width. 


An objective of Category 8 of the BANC-III workshop is to assess the solution variability caused by 
the propagation phase of CFD-based noise prediction. As such, a first step in assessing the variability is to 
compare the output from propagation codes provided the exact same input. Figures 14 and 15 show the 
difference between ANOPP2-F1A and upacs Acoustics predictions at the far- field and near-field microphones 
(using the impermeable surface as input) expressed as a ASPL about an average SPL of the two predictions; 
thus, the symmetric nature of the curves. Due to the high variance at extreme low and high frequencies, only 
frequencies between 0.24 kHz and 3 kHz were analyzed. Across the available frequencies and all observer 
positions, the variability of the predictions is less than 0.25 dB suggesting that there is good agreement 
between the two solvers when using impermeable data. Figures 16 and 17 show the far-held and near-held 
variability between ANOPP2-F1A, upacs Acoustics, and sAbrinA using the permeable data as input. Across 
all frequencies and observer positions, the variability is signihcant, particularly as the frequency increases. 
This variation in the high frequencies is partially due to the downsampled solution provided to sAbrina and 
suggests that the current analysis may be problematic because of the particular way the different propagation 
approaches were applied to this problem. 

While spectral deviation at all the observer positions as a function of frequency is interesting and allows 
an identihcation of key problem areas, it is advantageous to reduce this deviation to a single number, 
thus quantifying the variability in the propagation phase of noise prediction. The authors admit that the 
participation to Subcategory 2 is sparse and an estimation of the overall variability would benefit from more 
participation. Regardless, an overall estimate of the variability is calculated. Equation 11 shows the process 
for calculating the standard deviation across all microphones and frequencies (results from permeable dataset 
only). Taking the standard deviation for near- and far- field predictions results in a standard deviation of 1.5 
dB and 2.7 dB respectively. Due to the nature of the study, the authors feel that rounding up to 2 dB and 
3 dB for the near-field and far-held variability is appropriate at this time. This variability can be directly 
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Figure 13. Power spectral densities at microphone locations N1 through N5 and the origin from ANOPP2-F1A 
and upacsAcoustics codes using impermeable surfaces. Nine FFT averages were applied resulting in a 120 Hz 
bin width. 
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Figure 14. A SPL between ANOPP2-F1A and upacsAcoustics predictions at the far-field microphones using 
the impermeable surfaces as input. 
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Figure 15. A SPL between ANOPP2-F1A and upacsAcoustics predictions at the near-field microphones using 
the impermeable surfaces as input. 
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Figure 16. A SPL between ANOPP2-F1A, upacsAcoustics, and sAbrinA predictions at the far-field micro- 
phones using the permeable surfaces as input. 
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Figure 17. A SPL between ANOPP2-F1A, upacsAcoustics, and sAbrinA predictions at the near-field micro- 
phones using the permeable surfaces as input. 


attributed to the choice of propagation algorithm when performing a CFD-based noise prediction. 
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IV. Conclusion and Future Work 

This paper presents results from the BANC-III workshop, which explores variability in the propagation 
phase of CFD-based noise prediction. This included two test cases: an analytical solution of a quadrupole 
source near a sphere and a computational solution around a nose landing gear. Excellent agreement was 
found for the analytical test case between Integral Method solvers using the impermeable surface data. For 
the computational solution around a nose landing gear, differences of 2-3dB were common for impermeable 
surfaces. Potential causes of variability include error in the CFD solution provided as input, how the 
approaches used the CFD data, implementation details, and spurious signals in the IM solvers. Obviously, 
additional participation in future BANC workshops and additional categories will improve the assessment 
of the variability. 

This represents a first step in quantifying the importance of the propagation phase of CFD-based noise 
prediction. In the future, additional categories will be added, in particular ones that focus on calculating 
the scattering/shielding of noise and accounting for the influence of flow heterogeneities in the mid-field 
region. In addition to new categories, future work will involve modifying the resolution of the input data 
and adding more complex configurations of both Subcategory 1 and 2 that will allow a focused assessment 
of the variability due to influences such as computational reliability and flow characteristics. This will allow 
for improved understanding and assessment of the validity of predictions, including frequency ranges and 
observer positions. 
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